A BERNOULLI PROBLEM WITH NON CONSTANT GRADIENT BOUNDARY CONSTRAINT 
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ABSTRACT. We present in this paper a result about existence and convexity of solutions to a free boundary problem 
of Bernoulli type, with non constant gradient boundary constraint depending on the outer unit normal. In particular 
we prove that, in the convex case, the existence of a subsolution guarantees the existence of a classical solution, which 
is proved to be convex. 



1. Introduction 

' q ; Consider an annular condenser with a constant potential difference equal to one, such that one of the two 
i— plates is given and the other one has to be determined in such a way that the intensity of the electrostatic field 
^s constant on it. If O \ K represents the condenser, whose plates are CI and K (with K C CI), and u is the 
electrostatic potential, it holds Au = in Q \ K and |Du| = constant on either 3D or 3K, depending on which 
^rf them represents the unknown plate. 

This gives rise to the classical Bernoulli problems (interior and exterior), where the involved differential oper- 
j^ator is the Laplacian A, which expresses the linearity of the electrical conduction law. However, some physical 
situations can be better modeled by general power flow laws, then yielding to the p-Laplacian as governing 
^"operator. Moreover, one can consider the possibility to have a non constant prescribed intensity of the electric 

!^?ield on the free boundary. In particular, as the intensity of the electrostatic field E on an equipotential surface 

(Jj s related to its outer unit normal vector, through the curvature of that surface, one can assume | E | to depend 
+^)n the outer unit normal vector v(x) of the unknown boundary. In view of these considerations, we deal here 
rwith the following problem. 

i " Given a domain in CI C R N , a real number p > 1 and a smooth function g : S N_1 — > M such that 
^1.1) c < g(v) C for every V G S N_1 , 

^or some C > c > 0, find a function u and a domain K, contained in O, such that 

t> : 

<N ■ 
^1.2) 

on ■ 

o : 
o : 

^"j-tvhere v(x) = yk (x) is the outer unit normal to 9K at x 6 3K. 
> Here an later A p is the p-Laplace operator for p > 1, that is 

[Xj ■ ApU = div(|Du|P- 2 Du) . 

C$ If u is a solution to (|1.2|l we will tacitly continue u by 1 in K throughout the paper, so that a solution u to (|1.2| l 
is defined, and continuous, in the whole CI. 

The boundary condition |Du| = t has to be understood in a classical way: 

Jim |Du(y)| = |Du(x)|. 
yeo\K 

Moreover, in the convex case, that is when O is a convex set, we are allowed to consider classical solutions 
(justified by |H8l , since K inherits the convexity of CI, as shown later). 

Notice that, given K in (|1.2|l . the function u is uniquely determined since it represents the capacitary potential 
of O \ K; on the other hand, given the function u the free boundary 3K is determined as 3K = 3{x € M N : 
u(x) ^ 1}. Hence, we will speak of a solution to |[L2) referring indifferently to the sets K or to the corresponding 



Apu(x) =0 in O \ K, 

u = on 3D, 

u = 1, on 3K, 

L|Du(x)| = g(v(x)), on3K, 
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potential function u (or to both) and we will indicate the class of solutions as JP{C1, g), where O is the given 
domain and g is the gradient boundary datum. 

The original interior Bernoulli problem corresponds to the case p = 2, that is the Laplace operator, with 
constant gradient boundary constraint g(v(x)) = t. In general, given a domain CI C M N , and t > 0, the 
classical interior Bernoulli problem consists in finding a domain K, with KCQ and a function u such that 



(1.3) 



Apu(x) 







u 



u 




1, 



llDu| 



in D \ K, 
on 3D, 
on 3K, 
on 3K. 



Easy examples show that Problem (jl .3ft , and hence Problem (|1.2|l . need not have a solution for every given 
domain CI and for every positive constant T. Many authors consider the classical problem, both from the side 
of the existence and geometric properties of the solution. In particular we recall the pioneering work of Beurling 
[6 1 and several other contributions as HI, 0, |T2| , flT\. fT7\ . The treatment of the nonlinear case is more recent 
and mainly due to Henrot and Shahgholian (see for instance [14], [15]; see also [5],|8||, Ifl3l , 11201 and references 
therein). The uniqueness problem has been solved later in [9| for p = 2 and HI for p > 1. Here we summarize 
some of the known results. 



(1.4) 



Let O C I N be a convex C bounded domain. There exists a positive constant A p = A p (0), 
named Bernoulli constant, such that Problem (11.3ft has a solution if and only if t ^ A p ; in such 
a case there is at least one which is C 2 oc and convex. In particular for t = A p the solution is 
unique. 

In this paper we consider Problem (|1.2|l in the convex case, that is when the given domain is a convex set, and 
we prove that the convexity is inherited by the unknown domain without making additional assumptions on 
the function g. More precisely, let us indicate by ^~[C1, g(v)) the class of the so called subsolution to Problem 
(|1.2|l ; essentially, v and K are subsolutions if v solves 



A p v ^ 



inO\K 
on 3D 
on 3K ; 



= 

= 1, |Dv(x)| < g(v(x)) 

(see Section [23] for more details). 
Our main theorem is the following. 

Theorem 1.1. Let CI C R N be a convex C 1 domain, and g : S N_1 - 
IfJ?~{Cl, g[y)) is non empty, then there exists a C 1 convex domain K with KCQ such that the p-capacitary potential 
u of CI \ K is a classical solution to the interior Bernoulli problem (I1.2D . 



be a continuous function such that il.lt holds. 



The idea of a non constant boundary gradient condition has been developed in the literature by many au- 
thors, who considered the case of a space variable dependent constraint, a : CI — > (0, +00). We refer to 
|3|,[4j,[21j for a functional approach, and to 12, IZL EH for the subsolution method. In particular, an anal- 
ogous result to Theorem II. II has been proved in [ 16 J where the authors considered a Bernoulli problem with 
non constant gradient boundary datum a(x). For a given convex domain CI C M N , and a positive function 
a : CI (0, 00), such that 

c ^ a(x) ^ C, for every x£Q, 

for some C > c > 0, with 
(1.5) 

they consider the problem 



— convex in O, 
a 



(1.6) 



< 



f A p u(x) 
u = 
u = 1, 
[|Du(x) 



ax 



in O \ K, 
on 3D, 
on 3K, 
on 9K. 
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and they proved that, if a subsolution to the problem exists, then there exists a classical solution and moreover 
the convexity of the given domain transfers to the free boundary. 

Notice that in Problem dl.6ll the function a is required to be given in the whole CI, while in (|1.2|l the function 
g is defined only on the unit sphere S N_1 . Moreover, while in Problem (I1.6D the convexity property fll.51) is 
required for the boundary constraint a, in Problem (|1.2[) no additional assumptions on g are needed. 

2. Preliminaries 

2.1. Notations. In the N-dimensional Euclidean space, N ^ 2, we denote by | • | the Euclidean norm; for 
K C M N , we denote by K its closure and by 3K its boundary, while COnv(K) is its convex hull. Jf m indicates 
the m-dimensional Hausdorff measure. We denote by B(xo,r) the ball in M N of center xo and radius r > 0: 
B(xo, r) = {x G M N : |x — xd < r}; in particular B denotes the unit ball B(0, 1) and we set cum = Jff N (B). Let us 
define 

S N - X = 3B ={x G R N : |x| = 1}; 

hence ^(S^ 1 ) = Ntu N . 
We set 

m 

A m = {A = (Ai, A m ) | At ^ 0, Y_ A i = !}• 

i=i 

Given an open set CI C M N , and a function u of class C 2 {C1), Du = (u Xl , . . . ,u XN ) and D 2 u = (u XiXj )^ =1 
denote its gradient and its Hessian matrix respectively. 

2.2. Quasi-concave and Q 2 ^ functions. An upper semicontinuous function u : M N ->1U {±00} is said quasi- 
concave if it has convex superlevel sets, or, equivalently if 

u ((1 — A)xo + Axi) ^ min{u(xo),u(xi)}, 

for every A G [0, 1], and every xo, xi G M N . If u is defined only in a proper subset CI of M. n , we extend u as —00 
in M n \ CI and we say that u is quasi-concave in CI if such an extension is quasi-concave in M N . Obviously, if u 
is concave then it is quasi-concave. 

By definition a quasi-concave function determines a family of monotone decreasing convex sets; on the other 
hand, a continuous family of monotone decreasing convex sets, whose boundaries completely cover the first 
element, can be seen as the family of super-level sets of a quasi-concave function. 

We use a local strengthened version of quasi-concavity, which was introduced and studied in |19| : let u be 
a function defined in an open set CI C M. n ; we say that u is a Q 2 ^ function at a point x G CI (and we write 
ueQ^(x))if: 

1. u is of class C 2 in a neighborhood of x; 

2. its gradient does not vanish at x; 

3. the principal curvatures of {y G R n | u(y ) = u(x)} with respect to the normal — p^[ x j| are positive at x. 

In other words, a C 2 function u is Q 2 ! at a regular point x if its level set {x : u(x) = u(x)} is a regular 
convex surface (oriented according to — Du), whose Gauss curvature does not vanish in a neighborhood of x. 
By u G Q 2 ^ {CI) we mean u G Q 2 ^ (x) for every x G CI. 

2.3. Quasi concave envelope. If u is an upper semicontinuous function, we denote by u* its quasi-concave 
envelope. Roughly speaking, u* is the function whose superlevel sets are the closed convex hulls of the corre- 
sponding superlevel sets of u. It turns out that u* is also upper semicontinuos. 

Let us indicate by O(t) the superlevel set of u of value t, i.e. 

Cl(t) = {x G M N I u(x) ^ t}, 
and let CI* (t) = conv(Q(t)). Then u* is the function defined by its superlevel sets in the following way: 

0*(t) ={x£l N I u* (x) ^ t} for every t G M , 

that is 

u* (x) = sup{t G R I x G CI* (t)}. 
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Equivalently, as shown in IITOl , 




u*(x) = max |min{u(xi), ...,u(x N+1 )} : x t G CI \ K, 3A G A N+1 , x = 

Notice that u* is the smallest upper semicontinuous quasi-concave function greater than u, hence in partic- 
ular u* ^ u. Moreover, if u satisfies A p u = in a convex ring Q \ K (that is O, K convex with K C CI), then it 
holds ApU* ^ in CI \ K in the viscosity sense (see for instance |10|). 

2.4. Subsolutions. In his pioneering work (6J, Beurling introduced the notion of sub-solution for the classical 
Problem (I1.3D . This concept was further developed by Acker [1J and then generalized by Henrot and Shahgho- 
lian 1151 , |Tl6l to the case p > 1, both for constant and for non constant gradient boundary constraint. 

Following the same idea, let us introduce the class of sub-solutions to the generalized Bernoulli Problem 
H1.2)) . Let n be a subset of M N ; [CI, g) is the class of functions v that are Lipschitz continuous on CI and such 
that 

( A p v ^ in {v < 1} n O 

(2.1) < v = on 3D 

[|Dv(x)K g(v(x)) on3{v<l}nO. 

If v G [CI, g) we call it a subsolution. 

As in the definition of solutions, we say that a set K is a subsolution, and we possibly write K G J£"~ [CI, x) or 
(v, K) G &~ [CI, x), if K = {x G CI : v{x) ^ 1} for some ver {CI, x). 

In the standard case g = t, for some positive constant t, it is known that the class of subsolutions and that of 
solutions are equivalent, indeed in MT5H is proved that, if CI is a C 1 convex domain, and [CI, x) is not empty, 
then there exists a classical solution to 01.31 . In particular it is proved that 

K(Q,t) = C, u = sup v, 

ce^(o,T) vg^-(o,t] 

solve Problem dl.3l l and hence, recalling (|1.4|l , it follows as a trivial consequence: 

A p (Q) = inf{T : &{Cl,x) ^ 0} = inf{T : ^~{a,x) ^ 0}. 

Regarding the proof of Theorem ll.il it is clear that an analogous relation between subsolutions and solutions 
hold true also in the non constant case, that is: 

g)= |J C, u= sup v, 



solve Problem (|1.2|l and they are said maximal solution to (|1.2|) . 

3. Proof of the main result 

In order to give a proof of Theorem ll.il some preliminary steps are needed; they are collected in the following 
propositions and lemmas. 

Proposition 3.1. Let CI be a regular C 1 convex subset o/IR N ; let uo,ui G ^~{Cl, g) with Ko = {uo = 1} and 
K x = {ui = 1}. Define K = K U Ki, K* = convK. Then v G [CI, g), where v is the p-capacitary potential of CI \ K*. 
Moreover 

|Dv(x)| ^ g(v a(t) (u x )), 

for every x G CI \ K* and y x G dK* such that vk*(u x ) = — Dv(x)/|Dv(x)| =-v Q ( t )(x), being Cl[t) the superlevel set of 
v of level t = v(x). 

Proof. Let u* be the quasi-concave envelope of u = max{uo, ui}; it satisfies in the viscosity sense 

'A p u*>0 inO\K* 
= on 3D 

= 1 on 3K*, 
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and hence, by the viscosity comparison principle, 

(3.1) |Dv| < |Du*| on 3K*. 

Consider y G 3K*; then either y G 3K* n 3K or y G 3K* \ 3K. 

Assume y G 3K* n 3K, so that v K (y) = v K *(y). Then either y G 3K , or y G 3Ki and hence |Du*(y)| < 
|Duo (y )| or |Dui [y ) |; however in both the cases 

|Dv(y)| < |DiH(y)| < g(v K (y)) = g(v K .(y)), 

asu ,ui G ^~(Q,g). 

Now assume y G 3K* \ 3K. By Proposition 3.1 in [10] there exist Xi, ...,Xn G 3(Ko U Ki) such that Xi, ...,Xi G 
3Ko, %i+x, ...,Xn G 3Ki (with ^ I ^ N) and A G An such that 

VK (*i) = vx(xi) parallel to v Kl (xj) = v K (xj ) parallel to v K *(y), 

for i = 1, I, j = I, N and y = Xl£Li AiXt. Moreover thanks to Proposition 2.2 in (7| it holds 

where i k G {0, 1}. Hence, by d3Tb , ve^ (O, g). 

Notice that, as CI, K* are convex, the function v is quasi-concave, in particular, thanks Lewis's result [18], 
v £ Q-(^ \ K*). For every x G O \ K*, let Vm t ) (x) be the outer unit normal vector to the superlevel set 
My) ^ v(x)}; hence by Lemma 4.1 in HO, it holds 

|Dv(x)| < g(v K *(y x )), 

where y x G 3K* is such that vk» (y x ) = y o.{t) M- D 

For the sake of completeness we rewrite here two lemmas in [1611 which are particularly useful in the proof 
of Theorem ll.il 

Lemma 3.2 ([16]). Let Dr = {xi < 1} \ Br, where Br = B(xr, R) and xr = (— R, 0, . . . , 0). Assume 1 > and let ur 

sofoe 

ApU = in Dr 

u = I on {xi = 1} 

u = on 3Br. 

Then for any e > f//ere exists R sufficiently large such that |Dur| < I + e on 3Br. 

Lemma 3.3 ( |[T6l ). Lef u be the p-capacitary potential of the convex ring O \ K, with |Du| < C uniformly in CI \ K. Then 
any converging blow-up sequence 

u Tj (x) = — (l-u(TjX)) , 

af any boundary points gives a linear function uo = ooc[ , after suitable rotation and translation, where a = |Du(0)| 
and O indicates the origin. 

Following the idea of the proof of Theorem 1.2 in Ifl6l , now we present the proof of Theorem ll.il 

Proof of Theorem |LT1 Let us consider u = sup{v : v G [CI, g)}, and let u n be a maximizing sequence. Notice 
that, thanks to Proposition [3T] we can assume {u n } to be an increasing sequence of the p-capacitary potentials 
of convex rings CI \ K n , with |Du n (x)| < g("^K n M) on 3K n for every n. Let K be the increasing limit of K n ; 
hence K is convex and, as uniform limit of p-harmonic functions, u is the p-capacitary potential of O \ K, with 

|Du(x)| < g(v K (x)) on3K. 

We need to show that in fact |Du(x)| = g(vi<(x)) and we will prove it by contradiction, constructing a 
function w G J^"~ [CI, g) such that w ^ u with w > uat some point. Let us remind that v(x) indicates the outer 
unit normal vector to 3K at x. 

Let us assume by contradiction that there exists a point y G 3K such that 

oc= |Du(y)| < g(v(y)) 
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and assume y to be the origin with outer unit normal v parallel to the first axis. Let 5 be such that 

(3.2) <x + 3S<g(v). 
By Lemma 1331 the sequence 

U rj = — (l -u(TjX)) , 

converges to UoM = ood[ , hence for every rj > 0, 

(3.3) u(x) > 1 — otx^ — riTj , 

if rj is small enough, for x = (xi, ...,xm ) G B(0,rj). 
Consider 

Wr(x) =wr i£ (x) = <x + 



2) Va+6/2-e, 

where ur is as in Lemma [3721 and I = a + 5/2. Then there exist eo, Ro > such that for e < eo and R ^ Ro, 

(3.4) |Dw R | < a + 26, on 3{u R < e} = {w R = 0}. 

Moreover there exist 6i, 5 2 > such that 

w R >axj f + 8 2 on 3B(0, 1) n{xi > — 5i}, 

in particular we can fix Si small enough such that {u R = e} n 3B(0, 1) C {xi > —261}, and choose 

< 5 2 = 2inf{u R (x) - axi + : x G 3B(0,1) n{xi > -61}}. 

Let w(x) = 1 — rjWR(x/rj); notice that, as ur is quasi-convex, then w is quasi-concave. Moreover for rj suffi- 
ciently small, recalling d3 3| | it holds 

w < 1 — otx^ — 52Tj < u on3B(0,Tj). 

Define 



w x 



max{u(x), w(x)} inB(0,Tj), 
u(x) inM N \B(0,T i ), 



and W = {w = 1} = Tj{w R = 1}; observe that on 3B(0, rj ), w = w. By d372l and d374l , for every x £ Wit holds 

|Dw(x)| < a + 25 < g(v) -5. 
Notice that {ur = 0} = 3B(xr, R) and for every x G 3{ur = 0} it holds 

lim v Br (x) = v = (1,0, ...,0). 

R^oo 

Moreover lim E ^o{ u R = £} = { U R = 0} as limit in the Hausdorff metric of convex sets. Hence, by continuity, for 
sufficiently large R and sufficiently small e, we have 

|g(v) - g(v w (z))| 6, 

for every z G W n B (O, Tj ), and hence, 

|Dw(x)| < g(v) - 6 < g(v w (x)), 

for every x G W D B(o, rj). 

Then w G J£" _ (Q, g) and, since w > u at some points, we get a contradiction with the maximality of u. 
Therefore |Du| = g(v) on 3K. □ 
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4. Final remarks 

Remark 4.1. Notice that in the non constant case no characterization of functions g for which J^~(0, g) is 
not empty are known. However in some trivial case the existence or non-existence of a solution can be easily 
deduced by the characterization of the existence for the standard problem in (|1.4|) . Indeed if g satisfies 

min gM^Ap(O), then J*"(0, A p (0)) C J^(0, g), 

and hence &~ (O, g) ^ 0; on the other hand, if 

M = max i g(v) < A p (0), then &~ (O, g) C &~(C1, M) = 0, 

and hence problem (|1.2|l has no solutions. 

Remark 4.2 (Concavity property of Bernoulli Problems (O)). As in the classical case, geometric properties for 
the maximal solutions to (|1.2|l can be proved. Indeed, following the argument in (8j, it is possible to define 
a combination of the Bernoulli Problems (|1.2|l in the Minkowski sense and to prove that Problem (|1.2|) has a 
concave behaviour with respect to this combination. More precisely: fix A e [0, 1]; let Oo,Oi be two given 
convex domains and go, gi : S N_1 — > R + two continuous functions (which both stay far away from zero). We 
define Oa as the Minkowski combination of Oo, Oi, that is 

Q* = (1 - A)O + AOi ={ Z ={1- A)k + Aki | k e K , ki € K x }, 

and gA as the harmonic mean of go and gi, that is 

1 _ (1-A) | A 

9\{y) goM giM' 

Consider Problem (|1.2|l for Oo,go and Oi,gi, respectively; we define their combined problem of ratio A the 
Bernoulli problem of the type \\.2) . with given set Cl\ and gradient boundary constraint gA("v)- Following 
the proof of Proposition 7.1 in HQ we can prove that if JP~[Q.i, gt), i = 0, 1, are non empty sets, then so 
is ^~(OA,gA)- More precisely let (K(Oi, gi),ui) be the maximal solutions, for i = 0, 1 and let ua be the 
Minkowski combination of uo and ui of ratio A, that is 

{u A > t} = (1 - A){u ^ t} + A{m > t}; 

(see for instance |8[ for more detailed definitions and properties). The function ua belongs to JP~[Q.\, g\) and 
hence, by Theorem ll.il Problem jl.2} for Oa an d gA admits a solution (K(Oa, Qa^^a) which satisfies 

(1 - A)K(O , g (v)) + AK(0 1; 9l (v)) C K(O a , g A ). 

Remark 4.3 (A flop in the unbounded case). It could be natural to try to extend Theorem ll.ll to the unbounded 
case with an approximation method considering a sequence of given domains Or = Qn B(0, R) as R grows. 
As the sequence {Or} is monotone increasing by comparison principle K(Or, g) also increases and hence it 
converges to a convex set. Unfortunately, this approach fails in the limit process as it turns out that in fact Or 
converges to the given set O which means that the limit of maximal solutions degenerates. 

More precisely assume for simplicity O = M N , so that Or = Br = B(0, R) (or, analogously O = H~ the 
half space {xn < 0} and take Or = Br = B(xr, R), where xr = (0, 0, — R)). If R is sufficiently large, then the 
Bernoulli constant of Br, A p (Br) = Cn/R (see [8[ for example) is smaller than Co and hence 

B T CK(B R ,c ) CK(B R ,g(v)), 

where B T = B(0, r) is the unique solution to Problem (|1.2|l corresponding to O = Br and g(v) = A p (Br). 

Hence for sufficiently large R, ^~(Br, g) is not empty and Theorem 11.11 gives a sequence of quasi-concave 
p-capacitary potentials {u R } which solve Problem (|1.2)l in Or \ Kr, where Kr = B r . By easy computations 
one can check that r = R/cn, for some constant Cn depending on the dimension and hence, the sequence of 
interior domains {Kr}r>o is not bounded for R which tends to infinity. This implies that the limit of the maximal 
solutions (Kr, u r ) is not the solution to the limit problem. 

Acknowledgements. The author wishes to warmly thank Paolo Salani for his invaluable help, the several 
helpful discussions and, above all, for his encouragement and support. 
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